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Summary: Recent theoretical works on alternative metric theories of gravity give 
greater significance to solar-system tests of General Relativity. In particular, it is 
suggested that the post- Newtonian parameter 7 ought to be determined with great 
precision in order to discover possible preferred frame effects or effects related to a 
tensor-scalar theory of gravitation. In this work we focus on a future experiment, 
based on Doppler tracking of an interplanetary spacecraft near superior conjunc- 
tion. 

We simulate, under certain assumptions, the frequency of the signal received from 
a spacecraft equipped with the ultimate radio capabilities. The geometrical op- 
tics formalism adopted makes it particularly easy calculating the Doppler effect 
caused by a spherically symmetric refractive medium. This case includes both the 
post- Newtonian gravitational effect and a spherical solar corona model. Stochastic 
fluctuations in the plasma electron density, along with other sources of noise, were 
also considered in simulating the data. 

Applying a least squares fit to the simulated Doppler data, we are able to determine 
7 with a precision of 5 X 10""*, a factor four better than the present value. The 
availability of a two-frequencies link to the spacecraft and back would allow an 
improvement by at least an order of magnitude. Our results seem to indicate 
that measuring the time derivative of the Shapiro ranging curve is more promising 
than the ranging technique itself, although additional work comparing the various 
methods is needed. 



1. Introduction 

The Parametrized Post-Newtonian (PPN) formalism is the standard framework 
for performing experiments and formulating gravitational theories in the weak field 
and slow motion limit, appropriate to the solar system. In fact, in this limit, the 
spacetime metric predicted by nearly every metric theory of gravity has the same 
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structure. In other words, the metric can be expanded about the Minkowskian met- 
ric through dimensionless gravitational potentials of varying degrees of smallness. 
Metric theories differ one from another only in the numerical values of the coeffi- 
cients that appear in front of the metric potentials. The PPN formalism prescribes 
the values of these parameters for each particular theory under study. Usually, we 
consider only ten of these parameters, chosen in such a way that their values indi- 
cate general properties of the metric theories of gravity. Thus, the comparison of 
these theories with each other and with the experiments is straightforward in this 
framework [1]. 

After almost eighty years since General Relativity (GR) was born, we can con- 
clude that Einstein's theory has survived every test. This rare longevity, along with 
the absence of any adjustable parameters in the theory (namely two of the ten PPN 
parameters are equal to one, all others disappear), far from being satisfactory, seems 
to be the cause of a lasting search for a possible fatal discrepancy. Remarkably, as 
an upshot of these efforts most alternative theories have been put aside. Only those 
theories of gravity flexible enough to accommodate the resulting constraints have 
survived, the protective belt being provided by the free parameters and the coupling 
constants of the theory. 

The Brans-Dicke (BD) theory is certainly the most famous among the alterna- 
tive theories of gravity. It contains, besides the metric tensor g, a scalar field (f) and 
an arbitrary coupling constant lo, related to the PN parameter 7 by 

1 + ^ ... 
^=2T^- 

The present limit on 7 (see eq. (3) below) gives the very tight bound on the coupling 
constant u > 500. In addition, many Scalar- Tensor (ST) theories, generalizing the 
original BD theory, have been proposed in which the coupling function depends 
on the value of the scalar field, u> — u;{(f)). These theories are considered as the 
most interesting alternatives to GR, due to recent developments in cosmology (e.g. 
inflationary models) and elementary-particle physics (e.g. string theory and Kaluza- 
Klein theories). As recently pointed out by various authors (see [2,3]), a large class 
of ST theories contain an attractor mechanism toward GR in a cosmological sense. 
If this is what actually happened, then today we can expect, from eq. (1), a small 
deviation of 7 from the GR value 7=1. The present level of this deviation is 
estimated to be in the interval 

57 f« 10"^ - 10"^ (2) 



the exact value depending on the particular ST theory under investigation [3]. 
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2. A test based on Doppler technique 

The previous arguments provide strong motivations for experiments able to push 
beyond the present empirical accuracy on 7, which is [4] 

7 = 1.000 ±0.002, (3) 

where the quoted uncertainty is twice the standard deviation. 

A number of advanced experiments and space missions have been proposed, in 
order to improve this level of accuracy by at least two orders of magnitude (see, 
for example, [5]). In addition to these, the interplanetary mission CASSINI, which 
will be launched in October 1997, offers a new possibility for a test based on a free- 
flying spacecraft. The instrumentation will include a coherent Doppler capability 
with simultaneous tracking at X (~ SAGHz) and Ka (~ 33GHz) radio bands. Since 
CASSINI is the first deep space mission to implement Ka-band communications 
and tracking, a significant improvement over the best published results using S- 
band (~ 2.3GHz) and X-band is expected. 

In Doppler experiments, a radio link is transmitted from the Earth to the space- 
craft, coherently transponded and sent back to the Earth, where its frequency is 
measured with great accuracy, thanks to the availability of hydrogen masers at the 
stations used for interplanetary telecommunications. Comparing the transmitted 
and received frequencies, one can measure the Doppler shift 

Azv Idl 
uq c at 

where / is the overall optical distance (including diffraction effects) traversed by 
the photon in both directions. This technique has various applications, from the 
estimate of the gravity field of planetary systems (see [6,7] for recent examples), to 
the study of interplanetary plasma [8,9], and the attempts to detect low frequency 
gravitational waves [10-12]. Interplanetary Doppler tracking near solar conjunction, 
as a tool for testing relativistic effects, has been discussed previously in [13]. 

Besides the ordinary Doppler effect, of order v/c, two relativistic corrections of 
order (v/c)'^ contribute to the Doppler signal: the gravitational shift, which occurs 
when the gravitational potential at the source and the receiver is different, and the 
second order effect, well known in special relativity. Both effects have been tested 
using a space-borne hydrogen maser carried aboard a sub-orbital rocket, allowing a 
test of Einstein's prediction to about 10~^ [14]. A different, third order relativistic 
effect is observable in interplanetary Doppler measurements and, indeed, is implic- 
itly incorporated in the orbit determination numerical codes [15]. In experiments 
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near superior conjunction, with an impact parameter 6, the ray suffers a deflection 

4Mo 

° -8-10-^-^ , (5) 



where Mq = 1.5 km is the gravitational radius of the sun. This deflection changes 
the velocity component along the line of sight and will produce an additional 
Doppler effect of order 

v5=^v^8-10-''^. (6) 
b b 

Here and in the following we have taken G = c = 1. Because of the different 
technique, this can be considered as an additional test of GR, related to the time- 
delay effect in radar propagation flrst pointed out by Shapiro [16]. Since in a generic 
metric theory the expressions (5) and (6) are multiplied by {l+'j)/2, the experiment 
amounts to a measurement of 7. The current stability of hydrogen masers (10~^^ 
or better) would allow a signiflcant improvement in the accuracy of 7, assuming 
that the effect of the solar corona can be overcome using dual-frequency capability. 
Smaller effects, such as the gravito-magnetic force acting on photon propagation, 
can also be studied with a 'differential' Doppler technique [13]. 

Before the introduction of multi-link tracking signal, including links at Ka band, 
the technique of measuring the relativistic Doppler shift caused by solar gravity 
was not competitive with the ranging technique. Our initial assessment is that the 
Doppler technique will be more powerful than ranging for advanced experiments 
like CASSINI. The purpose of our work is to quantify that expectation. The plan 
of the paper is as follows. In sec. 3 we derive the gravitational Doppler effect, using 
standard geometrical optics formalism. The same formalism gives, in sec.4.1, the 
steady-state contribution from the solar plasma. Sec. 4. 2 describes how to simulate a 
real experiment. The simulated data are then analyzed in sec. 5, in order to predict 
the accuracy in the determination of the PN parameter 7, when single and dual- 
frequencies links are available. Some of the limits of our analysis are then presented 
and discussed. Finally, in sec. 6, we consider the Doppler technique in connection 
with two well-known alternatives, namely ranging and Differential VLBI (AVLBI). 



3. Gravitational Doppler shift 

It is well known that null geodesies of a static spacetime are curves of ex- 
tremal coordinate time interval between the two end-points [17]. In particular, null 
geodesies for the post-Newtonian metric 



-(1 - 2U)dt'^ + (1 + 2-^U)5^jdx'dx^ 
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generated by a massive central body (U = U (r) = Mq /r) can be derived from the 
variational principle 

6 J fidX = 0, 

where 

r 

is the refractive index, dX = 5ijdx^dx^ is the geometrical arc- length (in flat space), 
and the spatial path F is to be varied with fixed end-points. The photon's trajectory 
is thus obtained by integrating the ray-equation 

^(Aip) = V/i, (7) 

where p is the unit vector tangent to the path. Since we are dealing with the first PN 
effect, we can integrate eq. (7) along the zero-th order trajectory, i.e. the straight 
line from the source to the receiver. 



Fig. 1: The deflection geometry. 

We choose the geometry as in fig.l; 2; = is the plane through the sun, the 
earth and the spacecraft at any given time, and the unperturbed ray is directed 
along the x-axis, so that r = (A, b, 0). We get, in the limit £o,£i ^ b, 

+00 

A,-=yVMA=-2(ii^*. (8) 

— CO 

Projecting eq. (8) along the 7-axis we recognize the usual expression for the angular 
deflection of light, so that we will adopt the following helpful notation 

Ap = A9gr{b)y, Mgr{b) = -'^^^^-^^^. (9) 

The one-way Doppler shift is given by the invariant expression 

i/i _ vi • p(xi) 



vq Vo • p(xo) 



(10) 



where x^ and Vj are, respectively, the position and velocity of the receiver (emitter). 
To lowest order the photon momentum is p = (1, 1,0,0), which gives rise, when 
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substituted in eq. (10), to the ordinary special relativistic effect. This effect changes 
over the Keplerian time scale ^^P/Mq ~ i/t}. The gravitational contribution to 
the Doppler signal comes instead from the first order correction 6p to p, and it 
changes over the time scale b/v. These corrections are needed at the end-points. 
Since we are dealing with a static metric, the time component 5p^ = 0. The space 
components are obtained from eq. (8) after imposing the aiming condition, i.e. the 
passage through xq and xi: 

J SpdX^O. (11) 

-£o 

The result is 

Spi^o) = -Mgrib)j^y (12) 

IQ + ll 

Sp{^,) = +Ae,r{b)j^y. (13) 
Finally, substituting eqs .(12), (13) in eq. (10) we get the gravitational Doppler shift 

yAt) = -{ "'llf' )^0,r{b). (14) 

The first factor in eq. (14) depends on the spacecraft and Earth's orbits (according 
to fig.l, with we indicate the velocity orthogonal to the line of sight), while 
the second factor depends on time through the impact parameter b, which, near 
conjunction, changes almost linearly in time. Of course, the round trip, or two-way, 
Doppler signal is simply obtained multiplying eq. (14) by two. 

Figure 2 shows an example of this signal, as a function of time, when the 
spacecraft is on a circular orbit around the sun, with Mercury's orbital radius. We 
have assumed 7 = 1 in this case, and chosen the origin of time at closest approach 
{tc = 0). Moreover, the orbital plane is inclined by ~ 1° with respect to the ecliptic, 
so that the ray grazes the sun limb when b is minimum (in other words, b{tc) — Rq)- 
Note also that at closest approach the special relativistic effect, not shown here, has 
its minimum, since at this epoch the velocities are orthogonal to the photon's path, 
and we have only a second order effect. 



Fig. 2: The general relativistic Doppler shift due to the solar monopole. The inclination 

I between the ecliptic plane and the spacecraft orbital plane is large enough to prevent the 
solar occultation, i.e. I — 1^^ ^ b{tc) — Rq- The experiment is assumed to occur near 
the point where the satellite has the largest distance from the ecliptic. 
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4. Plasma effects 

The solar corona plasma is the main source of noise in a radio experiment near 
conjunction. In general, we can express the Doppler signal as 



y = ygr + Vpl + Vn , 



(15) 



where ypi is the dominant plasma noise, and i/n contains all non-dispersive sources 
of noise (clock, receiver, etc.), along with a minor dispersive contribution from the 
atmosphere. 

The plasma contribution ypi to the fractional frequency change (4) is related to 
the change in the optical path 



Al = 



where e is the electron's charge, me its mass, and A^g the total columnar content 
along the beam, A^e = J'T-e dX. Therefore, in order to calibrate for the plasma term, 
we should know the electron density along the path. We start by decomposing the 
electron density rig in a static, spherically symmetric part (ne)(r) plus a fluctuation 
Sriei i.e. 

ne{r,t) = {ne){r)-\-5ne{r,t). (16) 

As a matter of fact, the steady-state behavior is reasonably well known, and we can 
use one of the several plasma models found in the literature [18-21]. To be more 
explicit, we will refer to two particular models, namely-^ 



(ne)(r) = < 



2.99 1.55 
+ 



2.39 X 10^ 



,6 



X 10^ + 



3.44 X 10^ 



cm 



-3 



6 



+ 



1.67 X 10^ 



2.3 



cm 



model (a) 



model (b) 



(17) 



where 77 = r/RQ. The reason why we consider both models is two-fold: first, 
repeating the same analysis with two different electron distributions, we can test 
the model-dependency of the results presented here. The second reason is more 
subtle; in the following section we will describe a simulation and a fit of the real 
data, so that, if we want to obtain meaningful estimates for ^7, we cannot use the 
same model in both operations. More on this later. 



We do not consider here a correction factor which depends on the heliographic latitude. 
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4.1. Steady state coronal plasma 

We will now determine the Doppler effect due to the solar corona plasma for 
a model of the form (17a) or (17b). More precisely, we consider a generic plasma 
model of the form 

k 

According to the electromagnetic law of phase propagation through a plasma me- 
dium, the refraction index /j, can be expressed in terms of the plasma frequency i^p 
and the carrier frequency i/q as 



V V^^o/ 27Tmeiyi ^ ' 

We can now insert eq. (18) in the eikonal equation (7), and integrate along the 
unperturbed trajectory r = (A, 6,0), exactly as we did in the previous section for 
the relativistic ji (see fig.l). In the limit £0,^1 ^ ^ we easily find 

Ap = A^pi(6)y, 

where now 

and B{x,y) is the Euler's Beta function. Comparing A9pi with Adgr we notice the 
opposite sign - gravity bends the ray outwards, plasma inwards - and the different 
dependence on 6, plasma effect being steeper. The plasma deflection as a function 
of the solar offset b is shown in fig. 3. In this case, we adopted for the steady-state 
model the first of eqs. (17), and considered S, X, and K-bands radio frequencies. 
The absolute value of the GR bending is also shown. The same analysis for the 
second plasma model (eq. (17b)) does not yield any appreciable difference, so we 
shall omit it here. 



Fig. 3: The deflection angle as a function of the impact parameter b. The solid line is the 
plasma contribution, while the dashed line gives the absolute value of the general relativistic 
deflection, given by eq. (9) with 7 = 1. 

Imposing the boundary condition (11), as we did for the gravitational effect, we 
can now obtain 5p (x^). Inserting these in eq. (10) we eventually get the steady-state 



Relativity Experiments in the Solar System 



9 



part of the solar plasma Doppler shift. The details of this calculation have already 
been carried out in the previous section. The result is consequently analogous to 
the gravitational Doppler shift, i.e. 



Note from eqs. (14) and (20) that the first factor is the same as in ygr- Both effects 
are time dependent through the quantities £0,^1,^0,^1, and b. 

4.2. Towards a real experiment 

In the previous paragraph we have considered a spherical and static corona 
model, and derived the resulting Doppler shift, according to the law of geometrical 
optics. Unfortunately, the true electron density (16) contains also the fluctuations 
Sue, which require particular attention. In fact, these fluctuations are carried along 
with the solar wind speed V ~ 400 km/scc, so that their spatial scale is Vt ~ Rq/2, 
where r ~ 1000 sec is the temporal scale. On the scale h typical of the gravitational 
Doppler effect, one expects 5ne to be of the same order of magnitude as its average 



With the differential Doppler technique, introduced in [13], the plasma contribution 
can be greatly reduced by at least three orders of magnitude. However, this tech- 
nique requires a hydrogen maser aboard the spacecraft, which is beyond the scope 
of the present study. 

In conclusion, in order to be realistic we should consider a stochastic signal of 
the same order of magnitude as the the steady-state plasma effect ypi, as given 
by eq. (20). Therefore, we have generated a random gaussian noise of amplitude 
ypi{t), at each instant of time. For simplicity, we did not take into account the fact 
that this noise is not white; although this could be done quite easily, assuming for 
example a Kolmogorov's spectrum, we claim that this assumption does not affect 
very much our results. 

Finally, one must consider the non-plasma noise yn which contributes to the 
overall signal (15). To be conservative, we have assumed here 




(steady-state plasma) 



(20) 



[22], i.e. 



Sne{b,t) ^ {ne){b) . 



Vn 



~ 5 X 10 



-14 



(21) 



which takes into account the fact that the observations occur during day time, 
although the noise level (21) could be pessimistic in view, for example, of CASSINI 
specifications [23]. 
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Summing up all these contributions, we can give a reasonable description of the 
real data. These have been simulated in S, X, and K-bands, assuming the same 
geometrical configuration as in fig. 2, and adopting the solar corona model given by 
eq. (17a). The simulated data are shown in fig.4. A superficial look to the K-band 
plot could give an idea of the expected improvement over previous experiments, 
based on S and X-band tracking. Even a single K-band link could therefore provide 
a significant estimate of the PN parameter 7. The following section provides a 
tentative result for this expectation. 



Fig. 4-' The simulated data for the three bands. The axes and the flyby's geometry are the 
same as in fig. 2, where only the gravitational Doppler shift was shown. 



5. Regression analysis 

As we have seen, the main difference between the gravitational and the plasma 
Doppler effects is that the latter is much more steeper, as a function of b or t, than 
the former. This fact implies that we can try to isolate the gravitational contribution 
with a regression analysis. To be more explicit, given the N data points yi and their 
variance af, we want to minimize the function 



where j/det is the deterministic signal, sum of the gravitational term (14), with 7 
as a parameter, plus the steady-state term (20), which contains the parameters 
Q!i, Q!2, . . .. Thus 



We have implicitly assumed that the exponents (3k are not varied by the fitting 
process. As we anticipated, we have created the data yi using the three-component 
model (17a), but it would not be fair to use this information here, in the regression 
analysis. In other words, in order to take into account our ignorance about the 'true' 
spherical model that Nature has chosen, we will attempt a fit with the 'wrong' one, 
in this case the two-component model (17b), characterized by (3i = 6,^2 = 2.3. 
Despite its formal importance, this distinction is however irrelevant, as long as the 
models used are all reasonable approximations one of the other. 





(22) 
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We stress the fact that the orbital factor in the relativistic Doppler signal (14), 
is multiplied by a quantity of order 0{MQ/b), so that its Newtonian value suffices. 
The analogous quantity in the time delay expression (see eq. (26) below) is added to 
the relativistic correction, so that they are both required at the same order. This 
constitutes a major advantage of the Doppler method compared to the ranging 
technique. Nonetheless, it is obvious that the orbital quantities appearing in ydet 
are affected by statistical errors, due for example to non-gravitational forces. We 
will discuss these orbital stochastic effects, along with small variations resulting 
from dependence on the PN parameter in a subsequent paper. 

Before we proceed, we should address the problem of the optimal time sequence 
to be used. Since the plasma, both the steady-state and the stochastic terms, 
become more and more important approaching the closest approach (t = in fig.4), 
we can decide to cut-off that portion of the data, which introduces more noise 
than signal. Figure 5 represents the signal-to-noise ratio as a function of the cut- 
off impact parameter, indicated with bcut- Of course, dealing with the function 
avoids this complication. In fact, the post-fit residuals give ^7/7 as a function of the 
minimum impact parameter bcut considered (fig. 6). Comparing figures 5 and 6 we 
can conclude that the accuracy with which we can determine 7 increases as we track 
closer to the sun, up to a certain point (in the K-band case around bcut — S-R©), 
after which plasma effects begin to dominate and deteriorate the accuracy, so that 
we should not push the observations beyond this threshold. 

Fig. 5: The Signal- To-Noise Ratio as a function of the cut-off impact parameter bcut- For 
each value of bcut only the first portion of the data, from the beginning of the observations 
to the point where b — bcut> o-f^ considered. 

Fig. 6: The fractional error 5^/^ as a function of the cut-off impact parameter bcut > where 
bcut ho,s the same role as in Fig. 5. 

In conclusion, from our regression analysis we can infer that the availability of 
K-band tracking could allow us to determine the PN parameter 7 with a two-sigmas 
accuracy 

^ ~ 5 X 10-4 , (23) 
7 

a factor 4 better than the present value 2 x 10"^, but still far from the optimal target 
quoted in sec.l (see eq. (2)). However, this sensitivity would be enough to detect 
possible effects related to the existence of a preferred frame in the Universe. The 
solar system would have, with respect to this reference frame, a velocity ~ 10""^, 
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which therefore represents a natural threshold for possible effects at variance with 
GR. 

It seems very difficult to go much beyond the result (23) with a single link. 
Even in K-band, there are limitations imposed by plasma effects. However, since 
the plasma noise is dispersive, a drastic improvement can be achieved when simul- 
taneous multifrequency links are available. For instance, if we track the spacecraft 
in X and K bands, then from the two Doppler observables y^{t) and y'^{t) we can 
construct the quantity 

^"x' -^y,. + (24) 



which, in principle, does not contain the frequency-dependent plasma noise. A 
preliminary regression analysis shows that this observable offers the opportunity 
for a measurement of 7 with a precision at least an order of magnitude better than 
eq. (23). 

Given our simple simulation of the real data, these estimates represent our best 
guess about the outcome of a gravitational test based on Doppler tracking of a space- 
craft equipped with the ultimate radio capabilities. However, before we conclude, 
we would like to discuss the weak points in our analysis. To start with, we have 
simulated our data without imposing any strict obedience to an actual proposed 
experiment; although the opportunity offered by the CASSINI mission stimulated 
most of our work, we did not pay much attention to the CASSINI expected orbit, 
error budget, etc. [23]. When necessary, however, we always tried to assume real- 
istic, if not pessimistic, numbers. A more crucial caveat is that we neglected some 
possible complications, like the presence of the solar magnetic field, the evolution 
of plasma properties over the duration of the experiment, non gravitational accel- 
erations, and so on. In addition, our claim that we can get rid of plasma effects 
using multiple links is only partially true. For instance, one must consider that two 
radio signals at different frequencies traverse different paths, since the bending (19) 
is frequency-dependent. As a consequence, the plasma noise can not be completely 
removed from eq. (24). However, we can predict how relevant is this limitation. 
At closest approach, when the separation between the two frequencies is bigger, we 
have 

1/? 

A6o(fc) ^ ho{X, h) - h,{K, h) ^ , (25) 

where bo is the minimum distance of the ray from the sun (not to be confused with 
the impact parameter b) . Comparing this separation with the typical length scale of 
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the plasma clouds, Vr ~ Rq/2, we can conclude that this is not much of a problem. 
The result of this analysis is summarized in fig. 7, where the actual difference in the 
closest approach distance of the K-band and X-band radio signals is shown as a 
function of the impact parameter. 

Fig. 7: The difference Abo between the closest approach distances in X and K hands, as a 
function of the impact parameter b. The dotted line is the estimate (25). 



6. Future developments 

In the present work, we have considered a relativistic test based on Doppler 
technique. Of course, Doppler is not the only technique available when a free 

spacecraft is orbiting beyond the sun. Two other well established experiments may 
also be performed at conjunction, namely ranging experiments, aimed to observe 
the Shapiro time delay, and AVLBI measurements of the deflection of light by the 
solar mass. 

The Shapiro time delay [16] consists in the retardation of a light signal passing 
near a massive body (the sun in our case). The round trip light time (RTLT) from 
the Earth to the spacecraft and back is easily obtained from eq. (7) as 



Unfortunately, we do not have direct access to the Newtonian signal (At)^^^^ (this 
is somewhat related to the fact that we are using a coordinate-dependent quantity), 
so that what people do is a differential measurement of the variations in RTLT as 
the spacecraft passes through superior conjunction. This requires a good knowledge 
of the orbital motion of the target relative to the Earth. For this reason, the best 
result, quoted in eq. (3), has been obtained with the Viking landers, anchored to the 
surface of Mars [4], while the more recent free-flying Voyager 2 test gave ^7 ~ 0.03 



Improvements in the accuracy of Very-Long-Baseline-Interferometry (VLBI) 
made it possible to measure differential angular positions to the level of hundreds 
of microarcsec, making highly accurate measurements of the deflection of light pos- 
sible [25]. This experiment is conceptually simple, and does not involve secondary 
parameters which limit the confidence of ranging experiments. A series of observa- 
tions yielded a value |(1 -|- 7) = 1.000 ± 0.001 [26], comparable to the Viking test 
of the Shapiro time delay. 




(26) 



[24]. 
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In conclusion, every technique is based on different concepts, and presents its 
own problems and peculiarities. Therefore, deciding which one should be pursued 
during a given mission is a difficult task. The project's management has to make 
a decision based on the particular instrumentation available, on the designated 
orbit, etc. It might also be advantageous to combine the three techniques in an 
optimal way, keeping in mind the technical difficulties involved in a multiple-task 
experiment. Here, without attempting such approach, we have just considered one 
particular technique, and tried to define its possibilities. The comparison with the 
others will be the object of future studies. 
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